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Abstract 

We report on robust features of the longitudinal conductivity (Uxx) of the graphene zero-energy Landau level in 
presence of disorder and varying magnetic fields. By mixing an Anderson disorder potential with a low density 
of sublattice impurities, the transition from metallic to insulating states is theoretically explored as a function of 
Landau-level splitting, using highly efficient real-space methods to compute the Kubo conductivities (both axx and 
Hall axy)- As long as valley-degeneracy is maintained, the obtained critical conductivity axx — 1.4e^//i is robust 
upon disorder increase (by almost one order of magnitude) and magnetic fields ranging from about 2 to 200 Tesla. 
When the sublattice symmetry is broken, axx eventually vanishes at the Dirac point owing to localization effects, 
whereas the critical conductivities of pseudospin-split states (dictating the width of a axy = plateau) change to 
axx — e^/h, regardless of the splitting strength, superimposed disorder, or magnetic strength. These findings point 
towards the non dissipative nature of the quantum Hall effect in disordered graphene in presence of Landau level 
splitting. 

PACS numbers: 72.80.Vp, 73.63. -b, 73.22.Pr, 72.15.Lh, 61.48.Gh 



1 



Introduction. -The massless Dirac fermion nature of low-energy excitations in monolayer graphene 
remarkably manifests in the high magnetic field regime, where the energy spectrum splits up into non- 
equidistant Landau levels (LL) given by En = sgn{N) ^y2hvF'^eB\N\ |llil2|]. One fundamental signature 
of such peculiar spectrum is the existence of a fourfold degenerate zero-energy LL (twofold valley and 
spin degeneracies) where electrons and holes coexist. As a result, the integer quantum Hall effect (QHE) 
yD (measured in conventional two-dimensional electron gas) transforms to a half-integer (anomalous) 
QHE in graphene, with a quantized Hall conductivity given by cxxy = Ae^ /h x (A^ + 1/2) [4-^. 

Such anomalous QHE is tightly interwoven with the 7r-Berry phase and pseudospin degree of free- 
dom, and occurs as long as K and K' valleys remain decoupled [|7|]. In contrast, if disorder breaks sublat- 
tice symmetry and strongly mixes valleys, the QHE in disordered graphene is predicted not to differ from 
other two-dimensional systems, recovering axy = INe^ j h, with an integer [|7|-|9|]. Several experiments 
performed in high-mobility samples have revealed an additional quantized Hall plateau at a^y = 0, evi- 
dencing a splitting of the zero-energy LL which could result from spin and/or sublattice-degeneracy lift- 
ing, stemming respectively from Zeeman interaction, sublattice symmetry -breaking mechanisms [llOll . 



or electron-electron interactions QH 
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The origin of such quantized plateau at a^y = has been further discussed in relation with the 
measurement of a finite value of a^x at the Dirac point, suggesting an unconventional dissipative nature 
of the QHE, however difficult to decipher [Il5l-ll9ll. Indeed, a theoretical scenario proposes the existence 
of a dissipative QHE phenomenon near the Dirac point (with finite conductivity axx ~ 1 — 2eVh) which 
would be conveyed by counter propagating (gapless) edge states carrying opposite spin [Il5l-ll7ll. Finite 
axx (— e^/(7r/i)) at the Dirac point has been also obtained from numerical simulations in tight-binding 
models of disordered graphene (introducing either bond disorder [llSll or random magnetic flux [|20ll ). and 
related with the formation of extended states centered at zero energy (but in absence of a fully quantized 
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Differently, other experiments have reported a strongly divergent resistivity at the Dirac point 
(<^xx — )■ in the zero-temperature limit) [^)19il which has been analyzed in terms of a Kosterlitz-Thouless 
metal-insulator transition [|21|1 . The observation of a temperature-dependent activated behavior of axx (T) 
further points towards a non-dissipative nature of the plateau axy = for a spin- splitting gap opening 
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This shows that the experimental literature on QHE in graphene is very rich and diversified, mainly 
because there exist various qualities of graphene material (epitaxial, CVD-grown, or exfoliated from 
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graphite) as well as experimental measurement conditions (silicon oxide or boron-nitiride substrate and 
suspended graphene). Disorder can also exist in a large variety of flavors (adsorbed impurities, vacancies, 
grain boundaries) which thus demand for transport universalities to be established. 

In this Letter, the magnetic field dependent-fingerprints of the dissipative conductivity (a^x) of disor- 
dered graphene are explored with and without energy level splitting. Using a tight-binding Hamiltonian 
and real space order- quantum transport approaches, the Kubo conductivities axx and cr^y are com- 
puted as a function of disorder and magnetic field. By tuning the contribution of valley mixing, universal 
features of axx are unveiled, such as a robust critical conductivity of axx — 1.4e^//i at the Dirac point 
as long as valley degeneracy is unbroken. In contrast, if sublattice symmetry is lifted by some impurity 
potential, pseudospin-split states are generated and found to convey different critical bulk conductiv- 
ities axx — e'^/h, regardless the splitting strength and magnitude of the magnetic field. In between 
pseudospin-split critical states, axx is found to eventually vanish in the zero temperature limit owing to 
intervalley-induced localization effects, in conjunction with the appearance of the quantized axy = 
plateau. These findings establish different critical values of the dissipative conductivity at the center of 
Landau levels of lowest energies, together with a clarification on the non-dissipative nature of the QHE 
in disordered graphene in presence of energy level splitting induced by sublattice symmetry breaking. 

Methodology- Electronic and transport properties are investigated by using a simple tt-tt* orthogonal 
tight-binding (TB) model with nearest neighbor hopping 70 (taken as 2.7 eV) 

n = J2Vc.\a){a\ - 70 ^ e-"^-^\a){(3\, (1) 

a {a,l3) 

n 

where the magnetic field is introduced through a Peierls phase [[25] with a magnetic flux per hexagon 
being equal to (f) = <f A- dl = /^/e X] hexagon V^"/?- ^ suitable gauge is chosen, allowing the calculation of 
transport properties in disordered graphene with realistic values of B (here varied from about 2 to 200 
Tesla). An Anderson disorder is first introduced through a modulation of the potential profile, by taking 
onsite energies at random within [—W/2, W/2] (70-units) where W gives the disorder strength. This is a 
commonly used disorder model for explorin g th e metal-insulator transition in low dimensional systems 



(with or without applied magnetic field) [|29 



Quantum transport in high magnetic fields is studied with order-A^ computational schemes for 



axx{E, B) [|26|] . as well as for the Hall conductivity axy{E, B) 1281], using real space implementations 



of the Kubo approach. The scaling properties of axx can be followed through the dynamics of electronic 
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wavepackets using [|26h 



a{E,t) = e^p{E)j^X\E,t) (2) 

where p{E) is the density of states (DOS) and AX^(£',t) is the mean quadratic displacement of the 
wave packet at energy E and time t: 

Tr[(5(E - n)] 

A key quantity is the diffusion coefficient defined as Dx{E, t) = AX'^{Ef, t)/t, which gives the con- 
ductivity through Eq. (|2]) at a certain timescale. The spin degree of freedom is included as a factor of two 
for a and p, while calculations are performed with system sizes containing up to several tens of millions 
of carbon atoms and energy resolution down to 5 ■ 10"^7o. All the information about multiple scattering 
effects is contained in the time-dependence of Dx{E,t). The trace in Eq.© is evaluated numerically 
using random-phase wavepackets |v?rp) according to Tr [...] — )■ NMpRp\...\ipRv) Jldl- Such method has 
been now widely used for studying strongly disordered materials [|27|] . 

The Hall Kubo conductivity is also computed from the time evolution of random-phase wavepackets 
lyjRp) and the Lanczos method, by rewritting (Jxy{E) as 



n poo poo 

cTxy{E) = -- / dte-'^'"' / dE'f{E' - E) Re 

^ Jo J-oo 



{<^Rp\6{E' - H)]y— — \: — -3x{t)Ww>) 
E' — H + 17] 



, (4) 



with ix = ^[H, X], the current operator (X the position operator), while 77 — )■ is a small parameter 
required for achieving numerical convergence. A new algorithm has been implemented following prior 
studies [I28]. 

Critical conductivity a^xiE, B, W) of the zero energy LL.-We study the evolution of a^xiE, B,W < 
2.5) in presence of Anderson disorder (which preserves chiral symmetry). At the Dirac point, different 
transport regimes are identified in Fig|r|(main frame) (for t = 12 ps, the maximum computed time). 
At = 0, all states are localized by the magnetic field {axx tends numerically to zero), but small 
disorder brings delocalization into play, as manifested by the enhancement of axx with W. For non-zero 
disorder up to W=2, Dx{E,t) are found to saturate to some maximum values in the long-time limit 
(Dx{E,t) —7- D^a^(E)), pinpointing the establishment of a diffusive regime and absence of Anderson 
localization effects. Additionally, for small enough disorder (up to W c::^ 0.1), axx increases roughly 
linearly with W whatever the strength of the field (tuned from 4.8 to 100 Tesla). The value of daxx/dW 
depends on the magnetic field (being larger for lower B) as expected from the scaling of the magnetic 
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FIG. 1: (color online) Main frame: Zero-energy conductivity versus disorder strength W , and varying magnetic 
field (from 2.4 to 200 Tesla). Inset: Localization length for the zero-energy Landau level (solid line) and 
density of states (dashed line, and arbitrary unit) for 1^=1.5 and B=9.6 Tesla. 



length Ib oc B~^/'^, which suggests reduced disorder-induced delocalization effects as B is increased 
(and is shortened). 

A remarkable saturation of a^x to a constant value ^ lAe^ /h is further obtained for a large range 
of disorder strengths W E [0.3, 2] and magnetic fields varying between 2.4 Tesla and 200 Tesla (up to 
two orders of magnitudes). This value identifies the critical regime in which the interplay between dis- 
order and magnetic field preserves extended states only at the center of the LL, while remaining states 
become localized (a key ingredient of the QHE theory). This is further rationalized by analyzing the 
nature of electronic states in the vicinity of the Dirac point (for disorder W < 2). By converting the time 
propagation of wavepackets to their spatial spreading and infering a corresponding length-dependent 
conductivity axx{L), the localization lengths are extracted at selected energies by fitting cXxxiL) with 
an exponentially decaying function. The typical behavior of ^{E) is illustrated for W = 1.5 in Fig{T] 
(inset), with a diverging ^(E) with energy lowering. When disorder exceeds W=2.5, all states (includ- 
ing the states at the Dirac point) become localized and the system is driven to the insulating state (for 
experimentally accessible values of B) with disappearance of the QHE regime, in agreement with prior 
numerical studies yOfl. 

Critical conductivity of pseudospin split states. robustness of the obtained critical value at the 
Dirac point is further investigated by adding a density of impurities which break the local A/B sublattice 
symmetry. To induce pseudospin-splitting, we use a heuristic model which consists in shifting all onsite 
energies of A (and B) lattice sites by a constant quantity Va (and Vb). We first simplify to the situation 
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where all A and B sites are differenciated in energy according to Va = — Vb, which induces a splitting 
gap of Va — Vb = 2Va. The Anderson disorder potential is maintained but with | Va| ^ W, potentially 
masking the formation of a pseudospin-split gap (see uppermost curve in Fig. |2]inset). The superposition 
of both potentials mimics some weak imbalance in the adsorption site in the sense of a slightly preferred 
sublattice. We note that recent experiments curiously report such possibility of imbalance doping or 



structural damage [|31 
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FIG. 2: (color online) Density of states for various magnetic fields B and disorder strengths W and adding a 
superimposed A/B-sublattice impurity potential given by Va (see main text), which splits the zero-energy Landau 
level. Main frame: full coverage of sublattice impurities with pVa = Va =0.002. Inset: dilute random symmetry 
breaking potential with indicated percentage p (constant total strength pT^A =0.002) for two values of = 0.1 
(left) and W = 0.2 (right). 



Fig. [2l(main frame) shows the density of states for W = 0.2 and Va = 0.002 which corresponds to a 
weak imbalance of adsorption on one sublattice. This imbalance splits the zero energy LL for magnetic 
fields as low as B =9 Tesla, but the splitting is reduced with increasing field, and becomes hardly visible 
for 45 Tesla (close-up right part of inset, black). We note that the sequence of higher LL does not exhibit 
splitting for none of the studied magnetic fields. 

We next consider the situation where the imbalance potential between sublattices is diluted by adding 
Va only on a small percentage p of randomly selected A sites, while keeping the total strength pVa fixed 
for comparison (analogously —Va for randomly selected B sites with equal concentration). Note that the 
random uncorrelated part (characterized by W) always remains much stronger then the diluted correlated 
part. While for p = 100% every A site and every B site are shifted by Va and — Va, respectively, a 
lower value for p means a random distribution up to an extreme dilution of 0.1% (pVa = 0.002). The 



6 



corresponding DOS is displayed in Fig. |2] (inset) for a magnetic field of 45 Tesla and two values of W. 
It shows that even at low concentrations ofp= 1%, dilution has negligible effect on the DOS in terms 
of splitting and peak heights (a long as pV^ is kept constant). A further enhancement of the impurities 
dilution below p = 1% or increase of W (> 0.2) leads to a disappearance of peaks and splitting signature 
(Fig. [2l inset). Finally, we note that the splitting does not change the total weight of the DOS, i.e. the 
integrated DOS is unaffected by the splitting and the peak heights are half of the heights of the initial 
DOS. 

We then scrutinize the time dependence of a^x in the very dilute AB symmetry breaking potential and 
investigate how robust is the conductivity plateau seen in Fig. \T\ Fig. [3] gives the energy-dependence 
of aa:x at 25 Tesla for p = 0.1% (j)Va = 0.002, and W=0.2). One observes a broad feature for the 
conductivity at small times (t = 20), which does not show any zero-energy dip. This is consistent 
with the corresponding DOS (not shown) which also displays a single maximum (similar to the case 
of larger B in Fig. [2l inset). The axx at Dirac point is smaller than ^ lAe'^/h but more importantly 
displays a strong time dependence indicating the contribution of quantum interferences. Evaluating the 
quantum conductivity at short times (t = 0.07 ps) roughly corresponds to introducing an effective cut- 
off for quantum interferences, thus reducing localization effects. At longer times (t > 0.37 ps), this 
broad feature of the profile of axx is replaced by a double peaks structure, which stems from enhanced 
contribution of multiple scattering phenomena. Interestingly, the conductivities at the two peak positions 
(for electron-hole-symmetric points, indicated with dotted lines) are almost identical and marginally 
affected with time/length scales, which indicates that no localization effects develop at such energies, 
and corresponding (critical) states remain extended. 

Fig. |3](inset) shows axx{L) for several typical energies. For E2 (vertical dotted lines in Fig. |3](main 
frame)), axx{L,E2) ~ e'^/h and remain length-independent, locating the energy position of the new 
critical states at the center of the pseudospin-split levels. In contrast, axx{L, Ei) and crxx{L, E3) (peak 
tail) and axx{L, Eq) (band center) are seen to decay to zero, pinpointing the localization of corresponding 
states and transition to the insulating Anderson regime. 

The generality of our results is checked by performing a series of calculations for varying magnetic 
fields and different values for pVa which yields neither qualitatively nor quantitatively different results. 
Our main findings are summarized in FigSlat elapsed computational time t = 6ps and for W=0.2. The 
maximum value of the doubly peaked axx turns out to be 5-independent, which reminds the case of 
conserved AB symmetry (Fig. [U). In contrast however, two axx peaks are clearly observed. Surprisingly, 
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Energy (y^) 

FIG. 3: (color online) Main frame: Energy dependent conductivity for different times (p=0.1%, pVa =0.002, 
W =0.2). Dotted lines indicate positions E = ±pVa (which, for the case of p = 100% and W — s-O indicates the 
positions of split LL, cf. Fig. [U. Inset: axx versus length dependence for selected energies Eo=0, Ek=0.00l*k. 



the peak maxima a^x are not half of the maximum obtained in the unsplit case but reduced by a factor 
of ^ 0.7. This is a clear quantative difference which might be related to the massive/massless character 
of the Dirac electrons. Earlier works on the 2D electron gas have also debated on the critical value of 



dissipative conductivity p3L 134 ] . Fig|4] finally shows that (Jxx{E = 0) — while the double-peak height 
of ~ e^//i is robust for different magnetic fields and disorder strength pVp,- 

Finally, we scrutinize the evolution of the Hall conductivity axy at 45 Tesla for a weak and diluted 
potential (pVa =0.005, p=2.5%) that breaks A/B sublattice symmetry (inset of Fig. HJ. At the charge 
neutrality point, the zero-valued plateau axy = appears (black solid line) in contrast to the clean case 
(j)Va=0) where axy only crosses zero at a single point, when jumping from —2e^/h to +2e^//i (green 
solid line). The plateau width is here confirmed to be defined by the pseudospin-split states observed 
in the density of states. Note that in the case of ultraclean samples electron-electron-interaction effects 



have been found to also produce additional plateaus in cr^,^ [|14l1 . 

Conclusion - We have reported on robust transport features at the Dirac point for the zero-energy 
Landau level. In absence of energy level splitting, a critical conductivity axx — 1.4e^//i is obtained 
for magnetic fields ranging from about 2 to 200 Tesla. When A/B sublattice symmetry is broken by 
some imbalanced local impurity potential, pseudospin split states are found to convey different critical 
values axx — /h. A non-dissipative QHE is demonstrated in this model, since axx — > in between 
pseudospin-split critical states which further dictate the width of the axy = quantized plateau. Inter- 
estingly, very recent scanning tunneling microscopy experiments on intentionally chemically (nitrogen)- 
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FIG. 4: (color online) Longitudinal and Hall conductivities with and without sublattice impurity potential. Main 
frame: a^^{t=6ps, W = 0.2) using pVa=0.001 (dotted line), pyA=0.002 (dashed lines), and pyA=0.004 (solid 
Une). Magnetic fields as indicated. Left inset: (Ja;x(i=6ps) (solid line) and DOS (dashed line) for pVa = 0.004 
and W = 0.2. Right inset: Hall conductivity axy{W = 0) (for 6 = 0.0057o,r? = O.OOO570, soUd Unes) and 
corresponding DOS (dashed lines) for pVa=0 (green) and pVa=0.005, p=2.5% (black) at 45 Tesla. 



doped or hydrogen-functionalized disordered graphene have revealed the surprising manifestation of 



some sublattice symmetry 



of our findings plL 



32L 



35 



creaking mechanism, offering possibilities for the experimental confirmation 
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